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Real gas effects typical of the critical region of working fluids are a powerful tool to increase the energy 
performances of Stirling cycles, mainly at low top temperatures. To carry out the compression near the 
critical region the working fluids must have a critical temperature near environmental conditions and 
the use of organic working substances (pure or in suitable mixtures) as a matter of fact begins 
compulsory. The moderate thermal stability of the organic working fluids limits the maximum 
temperatures to 300-400 °C and as a consequence, the achievable cycles efficiencies result rather low. 
Carbon dioxide, with a critical temperature of 31 °C, is, among the traditionally inorganic gases, an 
exception and is considered here in comparison with organic substances. But the good thermodynamics 
of the cycles allows, in the considered cases, conversion efficiencies of about 20%, with good specific 
powers. The good energy performance of real gas Stirling cycles is obtained at the cost of high maximum 
cycle pressure, in the range of at least 100-300 bar. These high pressures nevertheless have large 
positive effects on the heat power transferred per unit of pumping mechanical power, and the low top 
temperatures have a positive influence on the material problems for the hottest engine parts. 

© 2010 Elsevier Ltd. All rights reserved. 


1. Introduction 

The Stirling engine works according with a closed and regen¬ 
erative thermodynamic cycle. In a traditional Stirling engine the 
volume variations are the consequence of the alternate motion of 
pistons carrying out the compression and the expansion of the 
working gas at two different temperatures. 

The cyclic variation of the volume produces mass fluxes inside 
the engine and the net effect is the production of useful work 
equals to the net heat exchanged, [1 j. 

The engine can use every kind of fuels. For example, natural gas, 
coal, biomass or solar energy, [2—5]. It is quieter than the traditional 
internal combustion engines. Owing to the internal regeneration, it 
may potentially reach very high thermodynamic efficiencies: the 
ideal cycle has in fact the same efficiency of the Carnot cycle 
operating between the same temperatures. In the ideal Stirling 
thermodynamic cycle the temperature increase of the working gas, 
passing from the cold region to the hot one, is induced by the 
regenerator. This is one of the main reasons that makes the Stirling 
particularly suitable for the realization of engines of small power 
(some kW, for example). On the contrary, in the thermodynamic 
cycles like the Carnot cycle, without regenerator, the increasing of 
temperature (and pressure) of the working fluid is totally assigned 
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to the turbomachineries, that are very difficult to realize if the 
engine is of very small size and power. 

Furthermore, the Stirling engines are in several aspects similar 
to the thermoacustic engines, [6], and, on the whole, the basic 
conclusions of the analyses carried out for a type are valid for the 
other too. 

The primary interest in the development of small power engines 
is in the distributed generation of heat and power, for example for 
domestic application. Small, cheap and highly reliable solar or 
biomass powered engines may even have a large market in the 
undeveloped countries. Often, an high specific power is a funda¬ 
mental prerequisite. 

In the field of the distributed power generation, as a thermal 
machine, the Stirling engine is in competition with Rankine cycles 
and micro-gas turbines. Among the advantages of the Rankine 
cycles is their good thermodynamic efficiency even at moderate 
high temperatures, provided that an adequate working fluid is 
selected; Brayton cycles, even in closed cycles, are also potentially 
very interesting but in any case the design of efficient small tur¬ 
bomachines results very difficult. 

In modern Stirling engines, the working fluids are usually 
helium or hydrogen for engines of high performances. Helium and 
hydrogen, by virtue of their good thermo-physical characteristics, 
allow the design of engines with high heat fluxes and moderate 
pressure losses, but on the other hand, they impose severe sealing 
problems. Air was repeatedly proposed for engines of simple design 
and relatively low specific power with low speed, [7,8] 
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Nomenclature 

<t> 

empirical parameter for the calculation of the pressure 




drops (Equation (8)) 



<P 

phase angle 

Symbols 

c P ,c v 

heat capacity at constant pressure or volume, J/kg K 

c0 

angular velocity (=27rN/60) 

m 

mass flow, kg/s 

Subscripts 

N 

revolution per minute, rpm 

cr 

critical point of the working fluid 

P 

pressure, bar 

max 

maximum 

Q 

thermal energy, J/cycle 

min 

minimum 

Q 

thermal power, W 

r 

reduced conditions (T r = T/T cr , p r = ptPcr) 

R 

gas constant, J/kg I< 

c 

compressor 

T 

(absolute) temperature, K 

E 

expansion 

t 

time, s 

R 

regenerator 

V 

volume, m 3 

H 

heater 

V swc , Vswe compression and expansion swept volumes 

I< 

cooler 

^cia V c i e compression and expansion clearence volumes 

T 

total 

z 

gas compressibility factor (=p/pRT) 

D 

dead 

w 

work, J/cycle 

CI< 

denotes a physical quantity on the interface between 

w 

mechanical power, W 


the compression volume and the regenerator 

w 

work parameter (=Wlp max V T ) 

1<R 

refer to a physical quantity on the interfave between 

£ 

fractional temperature difference (—min(e^, e H )) 


the cooler and the regenerator 

V 

cycle efficiency 

RH 

indicate a physical quantity on the interface between 

Vc, Ve 

compression and expansion polytropic efficiencies 


the regenerator and the heater 

V ideal 

ideal cycle efficiency (=1 -T K /T H = 1 - x) 

HE 

relative to a physical quantity on the interface between 

Vr 

thermal efficiency of the regenerator 


the heater and the expansion volume 

P 

density, kg/m 3 




As a rule, an increase of the mean operating pressure of the 
engine raises the useful work and an increase of the speed of 
the engine raises the useful power, but, if the engine is ideal and the 
working fluid is an ideal gas, the power parameter W* = Wlp max Vj 
is constant. 

Some Authors, [9], considered the possibility of using two- 
components and two-phases working fluids, aiming to a high 
specific output at moderate mean pressures. 

On principle, an increase of the specific work is also obtainable 
by means of chemically reactive working fluids, as, for example, 
dissociating gases, like N 2 O 4 , proposed in [10] for Brayton cycles. 

The Malone engine, [11,12], is the well-known regenerative 
engine, similar in all to the Stirling hot-air, but employing water as 
heat transfer medium instead of air. In the Malone engine the very 
low compressibility of the liquid led to very high maximum pres¬ 
sures (about 800—1000 bar), both the swept volume of the pistons 
and the rotational speed (30—300 rpm) were small. 

The potential advantages in using water as a working fluid are 
his good heat transfer coefficients and his high specific heat. The 
seal problems are likely easier to handle despite the very high 
pressures involved. 

According to [11,13], the furnaces of the Malone engines man¬ 
ufactured between 1925 and 1927 were stoked with coal and coke. 
Therefore it is likely that the temperature of the working fluid in the 
heater exceeded the critical temperature of the water (373 °C). 
Malone, [11 ], tested many liquids and, among them, carbon dioxide 
(critical temperature of about 31 °C). So, in any case, it is likely that 
the working fluid was a dense gas rather than a liquid. 

If an engine works with a dense gas instead of a liquid, at suit¬ 
able volumes variation and round per minutes, the maximum 
pressures are lower than that in the case of a liquid, and, at the 
same time, the power parameter W* = Wlp max V T would be high 
enough. On condition that a suitable working fluid is found. 

Moreover, if the engine operates in regions of the thermodynamic 
plane where severe real gas effects prevail, it is possible to obtain good 
efficiencies even at moderate high temperatures (with assured 


technological advantages), and where and when the working fluids 
behaving as ideal gases cannot assure acceptable performances. 

2. The numerical model and the assumptions of calculations 

There are many mechanical arrangements for the Stirling 
engines: the so called alpha—type, with two pistons; a piston and 
a displacer (in the same cylinder, the beta—type, or in two different 
cylinders, the gamma—type engine); some kinds of hybrid config¬ 
urations, as the Ringbom engine, [15]. Yet, in any case, the basic 
design parameters of the engine are: 

• the minimum operating temperature T R , at the cooler, and the 
hot temperature Th at the heater 

• the swept compression volume V svv c and the swept expansion 
volume V SW e 

• the total dead volume V D 

• the phase angle ip of the expansion space volume variations 
relative to the compression space volume variations 

• the speed of the engine N 

In the following, we will consider a two piston machine and 
a reference geometrical configuration with five different volumes, 
Fig. 1, [16]: 

1 . the compression volume Vc containing the working gas at the 
temperature Tc 

2. the cooler volume V R 

3. the regenerator, with a volume V R 

4. the heater volume V H 

5. the expansion volume Ve containing the gas at temperature Te 

The cooler, the regenerator and the heater are interposed 
between the two opposed pistons, coupled to a common crank¬ 
shaft, and the temperature of the gas, in the different operating 
volumes, is increasing from the minimum value (in the 
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Fig. 1. Mechanical configuration with five components and with opposed pistons used for the calculations. Adapted from [16], 


compression zone and in the cooler) to the maximum one (in the 
heater and in the expansion zone). 

The compression volume Vc and the expansion volume 14 
change in function of time t according to, [16] 


until the net thermal energy Qr exchanged by the gas in the 
regenerator on a complete cycle is zero. 

In the next subsections are illustrated the calculation 
assumptions. 


Vq — V dc + \V S wS 0 + cos(wt)] 

Vg = 14/e + jVsweP ® COS(wt + <p)] 

in which, w = 27rN/60 is the angular velocity; V dc and V de are 
respectively the clearance spaces in the compressor and in the 
expansion volumes; V swc and V SW e are on the contrary the two 
swept volumes. 

The total volume of the engine, Vj, is then assumed equal to 

V T = V clc + I4wc + V K + V R + V H + V de + V swe (2) 

and the dead volume V D is considered as 

Vd = V clc + V K + V R + V H + V de (3) 

The quantities to be evaluated at any instant of time, on a time 
interval equals to the period of the engine, are 


2.3. The thermal efficiency of the regenerator 

The regenerator is an essential component of the engine to 
obtain an high thermodynamic efficiency. Its effectiveness may be 
accounted by a thermal efficiency p R defined as the ratio between 
the heat transferred to the working gas over the cycle and the heat 
transferred to the working gas over the cycle when the regenerator 
is ideal, in ideal conditions the working fluid during his transit 
through the regenerator exchanges the maximum allowable 
thermal energy. In this situation one of the two temperature end 
differences, (Trk - TV) or {Th - Trh), is zero. 

If the regenerator is not ideal, there exists a minimum temper¬ 
ature end difference greater than zero and we can assume as an 
efficiency indicator, instead of p R , the design parameter 

e = min(£K,£ H ) (4) 


• the pressure p in all the compartments of the engine 

• the gas temperature 7c in the compression volume 

• the gas temperature T E in the expansion volume 

• the gas mass flow m C K from the compressor to the cooler 
(negative if the flow is from the cooler to the compressor) 

• the gas mass flow m KR from the cooler to the regenerator 
(negative if the flow is from the regenerator to the cooler) 

• the gas mass flow m RH from the regenerator to the heater 
(negative if the flow is from the heater to the regenerator) 

• the gas mass flow m HE from the heater to the expansion 
volume (assumed negative for the reversed flow) 


where 


e k 


Trk ~ Tk 
T h - T k 


and e H 


Th ~ Trh 
T h - T k 


(5) 


as a consequence, in the equations of conservation of the energy 
and of the mass, the temperature of the gas inside the regenerator 
and inside the cooler space have to be assumed respectively equals 
to T r = 0.5 x {T rk + T rh ) and to Tr= 0.5 x ( Trk + Tr ). The tempera¬ 
ture of the gas inside the heater is Th = 0.5 x ( Trh+T h ). 


2.2. Polytropic expansion and compression efficiencies 


Two further unknown parameters are the temperature T RK of 
the gas flowing from the regenerator to the cooler and the 
temperature T RH of the gas flowing from the regenerator to the 
heater over the cycle. If there are no pressure losses on the heat 
exchangers, the pressure p is the same in all the volumes; if the 
working gas is an ideal gas and the regenerator is ideal, then 
Trk = Tk and Trh = Th. 

The useful work over the cycle is W = W E - W c = 
f(W E -W c )dt and the cycle efficiency is evaluated as 
n = W/fQ H dt = W/Qh. 

The equations available are the usual equation of conservation 
of the mass and of the energy, applied to the whole machine and to 
every single volumes. The solution of the system of the differential 
equations gives the engine pressure p, the gas temperature Tc in the 
compression space, the gas temperature T e in the expansion 
volume and all the mass flows. Fixed the pressure drops and T rk , 
the solution of the system must be repeated at various Trh { <Th), 


The losses caused by fluid viscosity are accounted by means of 
the expansion and compression efficiency p £ and pc- In an ideal 
cycle r] E = ric= 1 . 0 ; in real processes of compression or expansion, 
according to the variation of the volume dV/dt, the expansion useful 
work and the compression necessary work result respectively 
lower and greater than the ideal ones. So, W E and W c begin 


... dV c 1 dV c n 

W c = p—,-S-— if —j-^<0 

c dt p c dt 

... dV c dV c 

W c = VcP 1[F >f ^>0 


( 6 ) 


• dV E dV E n 

We = VeP ^3F If ~dF > 0 

W £ = P ^1 if ^<o 

b dt p £ dt 


( 7 ) 
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2.3. The pressure drops 

For sake of simplicity, all the pressure drops on the heat 
exchangers are concentrated between compressor and cooler so 
that p K = p R = p H = p E and the mean pressure in the compression 
space pc, mean over the cycle results greater than the mean pressure 
in the rest of the engine. 

The pressure drops are first of all a function of the mass flows, 
but, for simplicity of calculation, we assumed the concentrated 
pressure drop constant over the time according to 

if m CK > 0.0 then p E = p c (\ — (j>N 2 ^ 

if rhcK < 0.0 then P£ = Pc(l+^N 2 j (8) 

in which <j> is a suitable empirical constant and N the speed of 
the engine. 

2.4. Further model assumptions 

The following further calculation hypotheses are assumed: 

• no leakage 

• compression and expansion are adiabatic (Q. c = Q £ = 0) 
processes 

• the working gas exchange thermal power with the surround¬ 
ings only in the cooler and in the heater 

• the external temperature walls of the cooler and of the heater 
are at the same temperature of the gas 

The variations of the gas density as a function of temperature 
and pressure and all the necessary thermodynamic properties are 
evaluated by means of the Peng—Robinson equation of state, [17], 
Chapter 8 , valid for pure fluids and fluids mixtures. 


v v-h v{v + b) + b{v - b) ^ } 

where 

b = £?=i 

a = £"=i £"= i x i x j \/afij ( 1 -%) 

b; = 0.07780 ^ 

Pcr,i 

a, = a,■(1)0.45724^ 

‘ «\ / Ptr.i 

<Xi(T) = [l +mi (l-yT ^)] 2 

m, = 0.37464 + 1.54226W; - 0.26992w? 

in which x; are the molar fractions, Iqj are suitable constants, and v, 
is the acentric factor of the component i of the mixture. For pure 
fluids, n = 1 and x, = xj = 1.0 and ky = 0.0. The parameter R is the gas 
constant. 

Then, from the volumetric behaviour, one can evaluate all the 
necessary thermodynamic properties of the working gas. 

3. Results of the thermodynamic analysis of the ideal engine 

The numerical procedure previously outlined was applied to 
helium (He), hydrogen (H 2 ), nitrogen (N 2 ) and argon (Ar), repre¬ 
sentative ideal gases at the ordinary temperatures, and to carbon 
dioxide (CO 2 ), trifluoromethane (HFC — 23, CHF 3 ), ethane (CHe), 
sulfur hexafluoride (SFe), pentafluoroethane (HFC — 125, CF 3 — 
CHF 2 ), difluoromethane (CHF 2 ). In Table 1 are some properties of 
the considered fluids. The parameters assumed for the calculations 


Table 1 

Molecular mass, critical temperature and critical pressure of the considered working 
fluids. 


Fluid 

Molecular 

mass 

Critical 

temperature 

Critical 

pressure 

(g/mole) 

(K) 

(bar) 

Helium 

4.003 

5.201 

2.275 

Hydrogen 

2.016 

33.25 

12.97 

Nitrogen 

28.01 

126.2 

33.98 

Argon 

39.95 

150.86 

48.98 

Trifluoromethane HFC-23 

70.01 

298.97 

48.36 

Carbon Dioxide 

44.01 

304.21 

73.83 

Ethane 

30.07 

305.32 

48.72 

Sulfur Hexafluoride 

146.06 

318.72 

37.60 

Pentafluoroethane HFC-125 

120.02 

339.17 

36.15 

Difluoromethane 

52.02 

351.26 

58.05 


are in Table 2 (a) and (b). An ideal cycle efficiency rather low 
(rudeai = 1 - t = 0.5) was assumed, to limit the maximum tempera¬ 
ture Th, on account of the thermal stability of the organic working 
fluids. 

The geometrical and kinematical parameters inTable 2 are just for 
reference and we have assumed those of Viebach CHP unit presented 
in [18], page. 213. The basic thermodynamic distinctive characteris¬ 
tics of the cycles depend in fact only on the thermal and volumetric 
comportment of the working gas. Even though, x, the geometric and 
the kinematic parameters control directly the value of W*. 

The critical temperature of the organic and inorganic considered 
fluids are such that they show severe real gas effects close to the 
standard environmental conditions. The critical pressure of CO 2 
may be perhaps too high for the particular application here studied. 
Anyway, it was considered as a typical case of study. 

Besides, carbon dioxide power cycles were extensively investi¬ 
gated in the past for applications where the final cooling of the 
working medium in a condenser takes place just below the critical 
temperature of the substance proposed, [22,23], and in recent years 
there was a revival in C0 2 cycles studies mainly connected with 
their perspective use in medium temperature advanced nuclear 
reactors, [24]. 

The refrigerant HFC — 23 was experimentally found stable in 
static tests at about 450 °C, [19] but its critical temperature of about 
26 °C implies a heat rejection process at relatively low tempera¬ 
tures, which could not be easy to perform in a hot environment. A 
second stable refrigerant, the refrigerant HFC — 125, with a critical 
temperature of about 66 °C is suitable for cogeneration systems. 


Table 2 

Assumed parameters for the calculations of the performances of the thermodynamic 
cycles. 


(a) 


^SWC = V swe 

0.5 dm 3 

Vclc = V c |e 

0.123 dm 3 

V K 

0.0944 dm 3 

Vr 

0.4054 dm 3 

Vh 

0.2934 dm 3 

Vd 

1.0392 dm 3 

Vj 

2.0392 dm 3 

r = T K /T H 

0.5 

N 

800 rpm 

<P 

7T/2 

(b) 


IDEAL CYCLE 


C3 

n 

II 

1.0 

e = min (e K ,e H ) 

0.0 

(j)N 2 

0.0 

REAL CYCLE 


C3 

n 

II 

rr? 

0.95 

e = min(£ K , e H ) 

0.05 

#1 2 

0.00055 
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The sulfur hexafluoride is a fluid widely used in industrial processes 
due to its high thermal stability and low toxicity, [20]. But it is 
a potent greenhouse gas and, associated with stainless steel, its 
thermal stability at high temperature is questionable, [21]. It is 
considered here because its critical temperature (about 46 °C) is 
intermediate between that of the ethane and that of HFC — 125. 

in any case, an effective way of obtaining a working fluid with an 
optimized critical temperature relying on the limited number of 
available substances is the use of mixtures, [14]. 

The Fig. 2.a represents the power parameter W* = Wlp max V T as 
a function of the maximum reduced pressure p max , r = PmaxlPcr for 
the working fluids listed in the Table 1. The cycles with ideal gas have 
W* = 0.035 at every p max , n the cycles working with gases displaying 
intense real gas effects at temperatures Tg have power parameters 
increasing from about 0.04 (in ideal gas conditions, at lowp maXi r ) to 
about 0.052 at p max , r ~ 1 1.5. At very high reduced maximum 

pressures (about 5.5 -h 6 ), the power parameter begins 0.1 h- 0.11: 
an increase, compared with that of the ideal gases, of 200 %. 

Calculated efficiencies for CO 2 ideal cycles are reported in the 
Fig. 2.b as a function of the maximum reduced pressure, for various 
reduced cooler temperatures Tk, r . As Tk. r increases the real gas effects 
become less intense and the cycle efficiency tends to the value of that 
of the ideal cycle using an ideal gas. For example, at p max , r = 2.0, the 
cycle efficiency is 0.26 for Tk, r = IT (power parameter W* = 0.058), 
and 0.44 for Tg, r =2.0 (power parameter W* = 0.04). Just for 
comparison, in the Fig. 2.b is reported the efficiency of hydrogen 
cycles at 'l' K =303.15 K. The efficiency is about 0.45 at every operating 
reduced pressure p max , r , and W = 0.035 (see the Fig. 2.a). 

The cycle efficiency in Fig. 2.b initially decreases with the 
maximum reduced pressure because real gas effects introduce 
a thermodynamic loss in the regenerator: although the temperature 
difference T RK - Tgis zero ( Eg= 0.0) the variations in the specific heat 
of the working fluid are such that Tfi - Trh results positive ( eh > 0.0). 



C 



For example, with reference of the curve at T K , r =Tg/T cr = 1.1 of the 
Fig. 2.b, the cycle efficiency is 0.26 at p m ax, r = 2, where the fractional 
temperature difference eh = {Th-T rh )I(Th-T k ) reaches its 
maximum value (see Fig. 2.c), and the cycle efficiency increases with 
Pmax, r because eh decreases. 

As a general rule, when the reduced cooler temperature 
increases the beneficial real gas effects on the net work become less 
important and the power parameter W* decreases. On the contrary, 
the thermodynamic efficiency improves. 

At very high reduced pressures, and at T K , r not too different from 
1 . 0 , the real gas effects are such that both the efficiency and the 
specific work are reasonably high. The maximum advantages are at Tg 
r near 1.0 and at high p max , r - In those conditions the power parameter 
and the thermodynamic efficiency actually may be very high. 

As an example, the cycle (1) in Fig. 2.d, at p max , r = 6.05 and 
Tg, r — 1-01, has an efficiency of 0.46 with e = ec = £h = 0 . 0 ; the cycle 
(2), at p max , r —1.06, has e c = 0 . 0 , but eh results 0.20 (,{Th-T rh ) = 
61 °C) and, as a consequence, a lower thermodynamic efficiency of 
0.23. In both cases, Th/Tk = 0.5. The power parameter is 0.1 for the 
cycle (1) and 0.05 for the cycle (2). In the Fig. 2.d is reported, just as 
a comparison, an ideal hydrogen cycle with Tg= 303.15 K. His 
efficiency is good, similar to that of the C0 2 cycle (1), but his power 
parameter W* results equal only to about 0.035. 

The remarkable difference in W* between real gas and ideal gas 
cycles is due to the high reciprocal isothermal compressibility 
(dp/dQ) r of dense gases at high pressures while the good efficiency 
of real gas cycles at very high pressure is a result of the distinctive 
variability of the heat capacity of the working fluid. 

In an ideal gas the heat capacities at constant volume C v and at 
constant pressure C p are only a function of the temperature while 
for a real gas they are also a strong function of the pressure. As 
a consequence, if the regenerator is ideal (as considered in our 
analysis) and the working fluid is an ideal gas, then Trk —> T K and 



d 



Fig. 2. Results for ideal cycles: t] C =t] E = 1.0, e = min(fjc, eh ) = 0.0, $ = 0.0.The cycles with carbon dioxide, ethane, HFC -125, HFC - 23 and sulfur hexafluoride are evaluated at Tk,t = 1.01; 
the cycles using nitrogen, hydrogen, argon and helium are calculated at T K = 303.15 K. (a) Power parameters W* as a function of the maximum reduced pressure for some working fluids, (b) 
Thermodynamic cycle efficiency as a function of the maximum reduced pressure at various T K , r for carbon dioxide and hydrogen, (c) Fractional temperature loss e H in the regenerator as 
a function of the maximum reduced pressure at various Tk, r for carbon dioxide, (d) Pressure-volume diagrams for carbon dioxide cycles at Tk, r = 1-01 and at two maximum reduced 
pressures and for hydrogen at p max = 100 bar. 












































































1708 


C.M. Invemizzi / Applied Thermal Engineering 30 (2010) 1703-1710 




Fig. 3. Cycle efficiency (Fig. 3 a) and power parameter (Fig. 3 b) for real cycles operating with various working fluids as a function of the maximum reduced pressure. The nitrogen 
cycle was calculated at Tk = 303.15 K. At this operating temperatures nitrogen behaves as an ideal gas. Cycles operating with the other fluids were calculated at Tk= 1.01 x T cr . The 
assumed polytropic efficiencies were t]c=Ve = 0.95. Pressure drops were accounted by the parameter 0 = 0.00055 IN 2 . 


Trh~>T h . If the regenerator is ideal but the working gas shows real 
gas effects, there are two possibilities: (a) the temperature Trk=Tk 
and, as a consequence, Trh < T H ; (b) the temperature T rk > T K and, 
as a consequence, T REI = Th- 

In other words, even if the regenerator is ideal, the compress¬ 
ibility effects on the specific heats produce a A T KR = Trr -T k > 0 or 
a ATjjh = Th- Trh > 0 in the regenerator: one of the two AT at the 
end of the regenerator is invariably greater than zero. 

However, at very high pressure the real gas effects are such that 
the heat capacities result chiefly, once again, a function only of the 
temperature, and as a consequence AT^^AT^ or, if the regen¬ 
erator is ideal, ATjjh = A Tkr = 0.0. 

4. Results of the thermodynamic analysis of the real engine 

The parameters r/c, r) E , <j> and e used to calculate the thermody¬ 
namic performances of real engines are reported in Table 2(b). 


Owing to the assumed simplifying pressure drops model, the 
pressure in the compressor space pc is greater than that in the 
cooler, regenerator, heater and in the expansion spaces. That is: 
Pc>Pk=Pr = Ph=Pe- 

In the Fig. 3 some calculation results of cycle efficiencies are 
reported for nitrogen (in Fig. 3.a), carbon dioxide, ethane, HFC — 23 
and HFC — 125. For the cycles operating with nitrogen, calculated 
at Tr= 303.15 K, the component inefficiencies are such as the cycle 
efficiency results practically zero at any maximum pressure (for 
example, 1.52% and 1.11% at pc, max, r = 6 and, respectively, e = ek = 0.0 
and 0.05). With the assumed parameters and in the considered 
range of reduced maximum pressure, the efficiencies of real cycles 
operating with helium and hydrogen result also practically equal to 
zero and are not reported in the Figure. Cycles operating with the 
other working fluids have been calculated at T Ki r = 1.01. 

In any case, at low maximum pressure pc, max the cycle effi¬ 
ciencies are very low (from 5 to 10%, depending on the fluid); 





DIMENSIONLESS TIME 


Fig. 4. Some calculation results for a cycle operating with ethane. The cooler temperature T K is 308.37 K (= 1.01 x T^), the heater temperature T H is T K /x = 616.75 K. (a) Time 
variations of the transferred mechanical powers, (b) Variations with time of the transferred thermal powers, (c) Dimensionless expansion and compression gas temperature as 
a function of the dimensionless time. 
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increasing the working pressure the cycle efficiency increases 
because the component inefficiencies become less important due 
to the very great positive compressibility effects on the thermo¬ 
dynamic of the cycle. 

The thermal loss in the regenerator, by itself, imposes a decrease 
in cycle efficiency of about 10—20 points, according to the fluid and 
the maximum operating pressure. A fractional temperature differ¬ 
ence e = e K of 0.05 corresponds to a minimum temperature differ¬ 
ence ( Trj(~Tk ) of about 15—17 °C, depending on the critical 
temperature of the considered working fluid. 

The Fig. 3.b represents the power parameter W* as a function of 
the maximum reduced pressure for four fluids working at Tg r = 1 . 01 . 
The optimum W* is 0.08—0.09 at about pc, max, r = Pc, max/Pcr=5 
for all fluids. 

The Fig. 4 refers to ethane, atpc, max. r = 4.29. The calculated cycle 
efficiency results 19% and the power parameter W* is 0.084. The 
mean pressures in the compressor space and in the expansion 
space are respectively 120 bar and 100 bar. The maximum ones are 
209 and 190 bar. 

The mechanical powers, normalized by the net power cycle, are 
reported in Fig. 4.a. The net power over the period changes between 
+9.3 and -5.8 times the net power cycle. The normalized thermal 
powers are in Fig. 4.b. The maximum regenerated power is about 70 
times the useful power, a further demonstration of the essential 
role of the regenerator also with real gas working fluids. Finally, in 
Fig. 4.c are reported the variations of the gas temperatures in the 
(adiabatic) compression and expansion volumes. The expansion 
temperature changes over the period of 22 %; the compression 
temperature changes of only about 11 %. 

5. Conclusions 


has a positive influence on material problems too: ordinary carbon 
steel could be an adequate choice even for the hottest engine part. 
However, the thermal stability could be negatively affected by 
contaminants, like lubrication oil. 

On the other hand, the good energy performance of real gas 
cycles is obtained at the cost of high maximum cycle pressure, in the 
range of at least 100—300 bar. These high pressures nevertheless 
have large positive effects on the heat power transferred per unit of 
pumping mechanical power. In fact, in regard to the heat transfer 
properties and to the power consumed by friction, we may observe 
that, as reported for example in [25], page. 445, the heat transfer 
coefficient h (the transferred heat power per unit of surface and for 
one degree of temperature difference) and the friction power Wf per 
unit of surface area may be represented by Equations (10) and (11). 

h = ^73^r r h(N Re ) ( 10 ) 

Np r 4r h 

(,1 > 

in which rj, represents the hydraulic radius, Ns e is the Reynolds 
number, Np r (= fiC p /k) is the Prandtl number, fz the gas viscosity 
and k the thermal conductibility. The functions F^Nge) and Fj(N Re ) 
are convenient functions depending of the geometry of the 
heat exchanger. 

Then, fixed the mean operating reduced pressures and 
temperatures, the geometry and the Reynolds numbers, the ratio of 
( h/Wf ) , for a gas which behaves as a real gas (compressibility 
factor z = z* less or greater than 1 . 0 ), and ( h/Wf) for the same gas 
in ideal conditions {z = z° = 1.0) results 


The results of the analysis carried out in the previous sections 
show that, above all, the real gas effects (in special mode the low 
gas compressibility at high reduced pressures) are responsible for 
very high power parameters W* in comparison with the values 
obtainable from cycles operating with ideal gases: 2—3 times, at T K 
near environment conditions, for example, for ideal cycle, at p ma x, r 
of about 5 and z = T K /T H = 0.5. 

On the other hand, for ideal cycles (r;c = Ve — 1-0, £ = </> = 0.0), the 
real gas effects introduce a thermal irreversibility in the regener¬ 
ator that lowers the thermodynamic efficiency. However, this 
negative effect decreases with increasing the maximum operating 
pressure. For example, the thermodynamic efficiency of the 
considered ideal cycles with carbon dioxide, at Tg, r = 1-1 is 0.25 at 
Pmax, r — P maxi Per — 2, and rises to 0.4 at p max , r = 6 (ideal cycle effi¬ 
ciency rn dea i = 1 - t = 0.5). 

So, at very high operating reduced pressures and with Tg, r = 1 
the real gas effects allow, in ideal cycles, high efficiency and at the 
same time high power parameters. 

The situation changes impressively introducing the component 
efficiencies: the real gas effects, although partly theoretically detri¬ 
mental in ideal cycles, give important efficiency gains in real cycles. 
For example, efficiencies of about 20%, with power parameters of 
0.07—0.09 are obtainable for cycles at Tg, r = 1-01, t= 0.5 and 
reasonable component efficiencies. At the same operating conditions 
(high i, relatively low regenerator effectiveness and relatively high 
pressure drops) engines with ideal gases, like hydrogen, helium or 
nitrogen, do not attain useful cycle performances. 

In order to make the critical region accessible to the compres¬ 
sion process, the working fluid critical temperature must have 
a well defined value: only organic fluids offer a wide selection of 
critical temperatures for different technical applications, but owing 
to the limited thermal stability of organic substances, only cycles at 
top temperature in the range of 300—400 °C are feasible and this 



( 12 ) 


for a fluid at mean reduced temperatures of about 1.5 and mean 
operating reduced pressures of about 5, the ratio of Equation 12 
may easily, roughly, be of the order of 10 . 

The design and manufacture of a machine operating at very high 
pressure (although at moderate maximum temperature) can pose 
technological problems, but the high W due to the low 
compressibility of the gas allows the realization of small engines 
with significant useful power and this (a) mitigates the structural 
problems, (b) permits the pressurization of the crankcase, with 
advantages in the seals realization. 
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